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Abstract 

The l/Nc rotational corrections to the axial vector constant and the 
isovector magnetic moment of the nucleon are studied in the Nambu - 
Jona-Lasinio model. We follow a semiclassical quantization procedure in 
terms of path integrals in which we can include perturbatively corrections 
in powers of angular velocity O ~ jy-. We find non-zero 1/Nc order correc- 
tions from both the valence and the Dirac sea quarks. These corrections 
are large enough to resolve the long-standing problem of a strong under- 
estimation of both QA and jJ^ in the leading order. The axial constant 
QA is well reproduced, whereas the isovector magnetic moment is still 
underestimated by 25 %. 
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Among the effective low-energy models of baryons the chiral soliton model, based on 
a semibosonized version IQ of the Nambu - Jona-Lasinio lagrangean 0, seems to play 
an essential role. In fact, an equivalent effective chiral quark-meson lagrangean can be 
derived 0] from the instanton model of the QCD vacuum. This chiral soliton model, 
which is frequently referred simply as NJL model, incorporates the general accepted 
phenomenological picture of the baryon as a bound state of A^^^ valence quarks coupled 
to the polarized Dirac sea {qq pairs). Operationally in the model the nucleon problem 
is solved in two steps In the first step, motivated by the large A*";, (number of 
colors) limit, a static locahzed solution (soliton) of a hedgehog structure is found. This 
hedgehog solution does not preserve the spin and isospin. Making use of the rotational 
zero modes in a semiclassical quantization scheme one can assign proper spin and isospin 
quantum numbers to the soliton. The model is rather successful in describing both the 
static nucleon properties and the nucleon form factors PJTO[]. The only exceptions 
are the axial coupling constant qa and the isovector magnetic moment /i^^, which are 
strongly underestimated {qa ~ 0.8 and /i^^ ~ 2.4(n.m.)) in the leading order, compared 
to the experimental values of 1.25 and 4.71 (n.m.). This is not a problem just of the 
NJL model. In fact, most of the chiral models, if they assume physical values of the 
pion decay constant /^r and the pion mass m^r, yield rather strong deviations from the 
experimental values for these quantities. In the case of the NJL model the simplest 
solution to this problem is to take into account the rotational corrections coming from 
the next to leading order terms. 



In a very recent paper Wakamatsu and Watabe [TT|] estimated in the NJL model 
1/Nc rotational corrections to qa thus providing a step in the right direction. They 
found a considerable non-zero valence contribution (^ 0.4). In this work they made 
the important observation that in the case of qa-, after the canonical quantization, the 
collective operators do not commute. However, in their scheme the non-zero result is 
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due to a particular order of the collective operators being not fully justified by path 
integrals or many-body techniques. As a consequence the valence contribution includes 
transitions between occupied levels, which violate the Pauli principle (even though, nu- 
merically, this Pauli violating contribution turns out to be only a tiny fraction of the 
valence contribution), and the Dirac sea l/N^ contribution vanishes exactly. The lat- 
ter is particularly puzzling since apart from the regular izat ion there is no any principle 
difference between the valence and the sea quarks in the NJL model. Based on the for- 
mulae of Wakamatsu and Watabe |]rT| Alkofer and Weigel studied the axial coupling 



constant in the context of the PCAC From their values for one can estimate how far 
the PCAC holds up in the linear order in VL and in particular, a violation of less than 2 
% for a reasonable value of the constituent mass M = 400 MeV can be guessed. 

Obviously a satisfying theoretical and numerical treatment of higher order rotational 
corrections in the semiclassical quantization scheme is still missing. It is therefore the 
objective of the present work to evaluate the 1/Nc rotational corrections to both the 
axial coupling constant and the isovector magnetic moment in the NJL model. To that 
end we will follow the theoretical scheme of Diakonov et al. elaborated in terms of 
path integrals. 

We start with the definitions of axial and isovector magnetic currents of a fermion 
field Al{x) = \I'+(x)7o7fc753^\l'(x) and Ji{x) = \I'+(x)7o7fc^^(a;), respectively. 

Here k means space components and a stands for the isospin index. We express the 
nucleon matrix element of a current operator as a path integral including quark 
\E', \E'"'" and meson fields U in Minkowski space: 



X 



(Ar(p/) |^+(o)6^^(o)| Ar(p)) = i / rf^xc/V"'"'' Vp-y 

T— »— ioo 

JVUJ V^V^+Jn{T/2, x)4(-T/2, y)^+(0)6^^(0)e^ / d^z^+D{U)^ 



The equality (|1]) should be understood as a limit at large euclidean time separation. Here 
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Z is the normalization factor which is related to the same path integral but without 
the current operator ^^O^^^ and 0% stands for the matrix part 7o7fc(75)^ of the 
current. The Dirac operator D{U) = idt — h{U) includes the single-particle hamiltonian 
h{U) = + M(3W^ + mo/3 with meson fields = e^^s"^ ' ^. Here a and (3 are the 
Dirac matrices and mo being the current quark mass. The composite operator: 

carries the quantum numbers JJ3, TT3 (spin, isospin) of the nucleon, where (3i is the 
color index, and ^ jj^ tt^ is a symmetric matrix in flavor and spin indices /j. 
In eq.(|l|) we can integrate the quarks out: 

(iv(p') |vi/+(o)d,-vi/(o)|iv(p)) = iriSiTT3riSlTT3A^c/ciWye-p'-e^p-y/pc/ 

X {(T/2,x|^|0,0)^,;,(6,"W(0,0|^| -T/2,y),,,,-Sp((0,0|^|0,0)O,") 

X (T/2,x|^| -T/2,y);,,J n(T/2,x|^| - T/2, y);,,.etr log/^(f/). (3) 

i=2 

In a natural way the result is split in a valence - the first term, and a Dirac sea contri- 
bution - the second one (see the diagrams on the l.h.s. of Fig.|l] a) and b)). 

In order to integrate over the meson fields U we start from a stationary meson 
configuration of hedgehog structure f7(x) = e^^ ' ^^(^) which minimizes the effective 
action. Then the integration over the meson fields U in the path integral can be done 
in a saddle point approximation, which is motivated by the large Nc limit. In the 
next step we should allow the system to fiuctuate around the static hedgehog solution 
U (x) making use of the rotational zero modes. Since the fiuctuations which correspond 
to the zero modes are not small they have to be treated "exactly" in the meaning of 
path integral. Operationally it can be done introducing a rotating meson fields of the 
form ?7(x, t) = i?(t)?7(x)i?^(t), where R{t) is a time- dependent rotation SU(2) matrix 

in the isospin space. It is easy to see that for such an ansatz one can transform the 
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effective action tilog D{U) = trlog(D(f/) — ^2) as well as the quark propagator in 
the background meson fields U 

(T/2,x|-^| -T/2,y) = RiT/2){T/2,^\^^Jy—^\ - T/2, y)i?+(-r/2), (4) 

where Q = —iR'^{t)R{t) = ^Qa'^a is the angular velocity matrix. Since f2 ~ (as can 
be seen below) one can consider as a perturbation and evaluate any observable as a 
perturbation series in Q which is actually an expansion in 

In this scheme the matrix element (eq.p of the current can be written as 

(A^(p/) |^+(0)6^$(0)| A^(p)) = ^NJ d^xd^yd^ze-'P'-''e'P ye'^P'-P>' JVR 
xD^^3^;(i?(T/2)){(r/2,x|^|0,z)i?+(0)6,»i?(0)(0,z|^| -T/2,y) 
- Sp((0, z 0, z)i?+(0)6,«i?(0))(r/2, X - r/2, y)} 

X (T/2,x|^| -r/2,y)2D^^3,3(i?(-r/2))etrM^-^). (5) 

Here, the finite rotation matrix D'^rp^j^, which carries the spin and isospin quantum 
numbers of the nucleon, appears due to the rotations R{t) of the valence quark propa- 
gators in eq. (|D correlated by the F jj^ matrices and the integral over z is due to the 
translational zero modes treated in the leading order. 

Now we are ready to make an expansion in Q. For the effective action, it is well-known 
procedure [|,^H,|lO|, which yields up to the second order in fl: 

6 r 

tT\og{D -Q) ^tilogD + i- J dtnl. (6) 

Here G is the moment of inertia. The first term will be absorbed in Z whereas the 
second one gives the evolution operator acting in the space of matrix R. Expanding the 
quark propagator 

1 1 1 o 1 

-Q h ... 
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we can separate the zero order (~ A'"^) and the hnear order (~ -^) corrections in Q. 
The expansion in Q is illustrated in Fig.[^ a) and b) for the valence contribution and for 
the Dirac sea one, respectively. 

Henceforward we will concentrate on the linear order terms. In this case we are left 
with the following path integral over R: 

B{T/2) ^ 

I VRDi^^j;iR{T/2))- Sp(i?+(0)r"i?(0)r'')r],(t)Z}:!^3,3(/?(-T/2))e^i 

R{-T/2) 

(7) 

Here we use the identity (i?+(0)6^i?(0))/g = | Sp(i?+(0)r"i?(0)r^)(6^)/g in order 
to separate the R{t) dependent part of the current which does not carry flavor indices 
fg. The path integral (^) can be taken rigorously within the approximation (|^). We 
obtain the well-known canonical quantization rule Qc Jc/Q, where Ja is the spin 
operator. The final result for path integral is a time ordered product: 

^{-t)Dab{RmJc + mJcDabiRm. (s) 

which should be sandwiched between the nucleon rotational wave function. In order to 
obtain the result we essentially made use of the basic feature of the path integral: 

J VqFMti)) ■ ■ ■ Fn{q{tnW' = (g2, ^ |r{A(g(ti)) ■ ■ ■ Fn{q{tn)}\ gi, Ti), 
qi=giTi) 

namely that the path integral can be equivalently written as the expectation value of 
the time ordered product of the corresponding operators. Using (|^) and the standard 
spectral representation of the quark propagator it is straightforward to evaluate the 
matrix element of the current eq.(|^). In particular for the linear correction we get 

{N{p') \^^0)dt^m Nip)f' = (J, J3T3 |[^, D,,]\ J, J3T3) 
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xN^Y. -^—{m\Ta\n) /rf=^^e^(p'-P)-$+(z)6^$„(z). (9) 

n>val 
m<val 

Here and are the eigenfunctions and eigenvalues of the single-particle hamiltonian 
h. Since the r-matrix element is asymmetric with respect to exchange of the states m 
and n a commutator appears in the collective matrix element of eq.(|^). It can be easily 
calculated: 

(J = 1/2, ^3^3 I Dab]\ J = 1/2, J3r3) = -^^£e635a3. (10) 

Using that both the axial coupling constant and the isovector magnetic moment are 
related to the corresponding form factors at = finally we get for the 1/Nc rotational 
corrections: 

9a{^^) = ^7^ E -^—{m\Ta\n){n\[(TXT]a\m) (11) 

^ n>val ^" ^"^ 

m<val 

and 

/^(^^) = ^l^ E {m\Ta\n){n\j5[[(TXx]xT]a\n) (12) 

^ n>val ^™ 

m<val 

In both eqns.(|ll|) and (0), a summation over a is assumed. As should be expected both 
expressions have similar structure with transitions from occupied to non-occupied levels 
and back. They include also an essential non-zero contribution from the Dirac sea. In 
contrast to the leading order the above expressions are finite and one does not need to 
regularize them. 

The parameters of the model are fixed in the meson sector to reproduce = 93 
MeV and m.^ = 139.6 MeV. Similar to other works @,PJ^ we use a numerical procedure 
based on the method of Kahana and Ripka ||14[. The results for qa and fi^^ up to the 



first order terms in Q, are presented in Fig.^ and Fig.^ as a function of the constituent 

quark mass M. As can be seen in leading order (Q^) they are almost independent of 
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the constituent quark mass and the valence contribution is dominant. In the next to 
leading order (Q^) the valence and Dirac sea contributions show much stronger and quite 
different mass dependence: with increasing constituent mass M the valence part gets 
reduced whereas the Dirac sea part increases and becomes dominant. However, their 
sum shows almost no dependence on M. The result of Wakamatsu and Watabe for qa 
(labeled as W&W in Fig.|^) deviates from the present numbers. This is due to the fact 
that in their scheme the contribution of the Dirac sea to qa in the linear order in Q 
vanishes exactly. In the present calculations for both quantities the enhancement due to 
the l/Nf. rotational corrections improves considerably the agreement with experiment. 
In the case of qa the experimental value is almost exactly reproduced and the inclusion 
of next order corrections will perhaps even overestimate it. In contrast to qa for the 
isovector magnetic moment /i^^ we are still below the experimental value by 25 %. It is 
interesting to notice that the enhancement for both quantities due to the 1/Nc rotational 
corrections is very close to the estimate ^^f^ ■ 

To conclude, we have evaluated the axial coupling constant qa and the isovector 
magnetic moment fi^^ in the Nambu - Jona-Lasinio model in the next to leading order 
in the semiclassical quantization scheme. The 1/Nc rotational corrections are large 
enough to resolve the problem of strong underestimation of these two quantities in the 
leading order. In particular, qa is almost exactly reproduced. However, such large linear 
order corrections imply that in order to control the perturbation series in Q the next 
order (l/N^) corrections should be investigated as well. 
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FIGURES 



FIG. 1. Diagrams corresponding to the expansion in Q of the current matrix element: a) 
the valence contribution and b) the Dirac sea contribution. 

FIG. 2. Axial vector coupling constant gA evaluated up to the linear order in O as a function 
of the constituent quark mass M. 

FIG. 3. Isovector nucleon magnetic moment evaluated up to the linear order in CI as 
a function of the constituent quark mass M. 
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